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Abstract We introduce a potential of multipartite entanglement for a 
system of n qubits, as the average over all balanced bipartitions of a bipartite 
entanglement measure, the purity. We study in detail its expression and look 
for its minimizers, the maximally multipartite entangled states. They have 
a bipartite entanglement that does not depend on the bipartition and is 
maximal for all possible bipartitions. We investigate their structure and 
consider several examples for small n. 
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1 Introduction 

Entanglement is one of the most striking features of quantum phenomena 
[¥T] , It plays very important roles in quantum information processing such 
as quantum computation [37], quantum teleportation [B] (for discussions 
on experimental realizations see [11X12, 25. 38]). dense coding [JJ and quan- 
tum cryptographic schemes (TTI[T5ll2"3] . Nevertheless, the quantification of 
multipartite entanglement is no simple matter. 

Entanglement is intimately related to the very mathematical structure 
of quantum mechanics and complex Hilbert spaces. In particular it is a 
straightforward consequence of linearity (superposition principle) in tensor 
product Hilbert spaces (composite quantum systems). 

Consider a quantum system composed of two parts (e.g. two particles): 
part A, whose Hilbert space is Ha, and part B, whose Hilbert space is Hb- 
According to quantum mechanics, the composite system lives in the tensor 
product Hilbert space H = Ha ® Hb ■ The most familiar example is that 
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of two spinless particles, whose Hilbert space is L 2 (R 3 ) eg) L 2 (R 3 ) = L 2 (R 6 ). 
The linearity of TC implies that the states \ip) of the composite system Tt 
are linear combinations of product states, namely, 

W) = ^ZijlVi) ® \Xj), 

ij 

with \ifi) G TLa and \xj) € T~Cb- This entails interference among probability 
amplitudes of two-particle states, that is the analogous of the simpler case of 
one-particle interference. For example, the probability amplitude of having 
both particle A in state \ipi) and particle B in state interferes with the 
probability amplitude of having both particle A in state \(f2) and particle 
B in state |x2)- As a consequence there exist correlations of quantum na- 
ture -entanglement- between quantum subsystems. These correlations are 
stronger than the classical ones, in the sense that they violate a class of 
inequalities, named after Bell, that must be satisfied by all classical corre- 
lations [31] . 

The most striking violation of Bell's inequalities is given by a particular 
class of states: maximally entangled states. The simplest example is that of 
two spin- 1/2 systems (or qubits), whose Hilbert space is C 2 Cg> C 2 , in the 
singlet state _ 
|#) = (|0)®|1)-|1>®|0»/V2, 

where {|0), |1)} is the natural basis of C 2 , representing spin up or down in 
a given direction. The expectation value of any local observable O of the 
first spin is given by 

<*|0 ® 1\<P) = l(0\O\0) + ~<1|0|1) = Tr(p A 0), 

and thus is an incoherent average corresponding to a completely mixed 
reduced density matrix of the first spin pa = 1/2. Analogously for the second 
spin. Therefore, spin measurements in a given, arbitrary, direction over an 
ensemble of pairs prepared in a singlet state will result in a completely 
random sequence of and 1. On the other hand, the results of joined local 
measurements exhibit strong correlations, due to the fact that the total spin 
is 0: the two spins are always found pointing in opposite directions. The two 
random sequences arc exactly complementary. 

Maximally entangled states are characterized by the property, just shown 
for the two-qubit singlet state, that to a perfect knowledge of the state of the 
composite system corresponds a complete ignorance of the states of its two 
parts. More precisely, although the composite system is in a well determined 
pure state, its two parts are in completely mixed states. See Corollary [3[ 
Therefore, all information is totally shared by the two parts. Note that this 
situation is strongly at variance with the classical case, in which a complete 
knowledge of the total system is equivalent to a complete knowledge of both 
its parts. In quantum mechanics this is only a necessary condition. 

In general, the degree of bipartite entanglement of a quantum system 
can range from a maximum, when its two parts are in completely mixed 
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states, to a minimum, when its two parts are in pure states, and thus a 
complete knowledge of both subsystems implies a complete knowledge of 
the entire system, as in the classical case. This is the case of separable, or 
unentangled, states = \ip) ® |x), that have no correlations between the 
two parts and, thus, no shared information. 

The degree of bipartite entanglement [48] of a composed quantum system 
can be consistently quantified, as in Definition [5J in terms of the purity of 
the reduced density matrix of one of the two subsystems (purity can be 
proven to be the same for both, see Lemma [5]). A lower value of purity will 
correspond to a larger value of entanglement. 

On the other hand, there is no unique way of quantifying multipar- 
tite entanglement [4] , that is entanglement among n given parties of a given 
quantum system. Different definitions often do not agree with each other, be- 
cause they adopt different strategies, focus on different aspects and capture 
different features of this quantum phenomenon [14, 15,33,36,47]. There is a 
profound reason behind this manifest disadvantage: the number of real num- 
bers, i.e. the invariants under local unitary transformations l,2,3j, needed 
to quantify multipartite entanglement grows exponentially with the size of 
the system, so that the definition of appropriate entanglement measures, 
able to summarize the most salient global features of entanglement, can be 
very difficult. 

A natural generalization of the bipartite case is to quantify the entan- 
glement among n parties by considering the average purity over subsystems 
[211I22| . In this paper we will follow this strategy and, in particular, we will 
consider systems of n qubits. 

After introducing notation and discussing some results about bipartite 
entanglement in Section[2j we move to multipartite entanglement and study 
the properties of the potential of multipartite entanglement (i.e. the average 
purity over balanced bipartitions) , Definition 1 101 and of its minimizers, i.e. 
quantum states with the maximal degree of multipartite entanglement. In 
the ideal situation, the bipartite entanglement of such states is not only 
maximal, but also does not depend on the way one decides to bipartite the 
total system into two subsystems. See Definition [9l Our approach is based 
on the action of the permutation group on the Fourier coefficients of the 
quantum state and thus is of combinatoric nature. 

The potential of multipartite entanglement, Eq. ([32)) of Theorem [9[ is a 
quartic Hamiltonian 



where z — (zk), with k = (fci, • • • , k n ) E {0, l} n is the vector of the Fourier 
coefficients of the state of a system of n qubit in the computational basis 




k,k',l,l' 
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One of the aims of this paper is to investigate the complex structure of the 
long-range coupling function A, that appears in the above expression. This 
is accomplished in Theorems [T0j[TT] and [13l and in CorollaryO In particular, 
a measure of the complexity of the potential of multipartite entanglement 
is given by the number of its nontrivial interfering terms, that scales like 
2"~ 3 3™ (see Theorem [T2l and the following remark). 

Sections HI and [5] will then be devoted to investigate maximally multi- 
partite entangled states (MMES) [31], i.e. the minimizers of the potential 
of multipartite entanglement. In particular, the structure of perfect MMES, 
i.e. minimizers that are maximally entangled with respect to any biparti- 
tion, is analyzed in Section [4] by making use of a probabilistic approach, 
Theorem [15] gives a complete characterization of their population probabil- 
ity vectors (|zfe| 2 ), while Theorem [TBI exhibits the equations that must be 
satisfied by their phases = {zk/\zk\)- The number of equations quickly 
overcomes the number of variables, since their ratio scales as 2 n+1 /y / n with 
the number n of qubits. See Theorem 1171 and following remark. Therefore, 
for large systems it becomes more and more difficult to have a perfect MMES 
solution, unless symmetries subtly conspire to reduce the number of inde- 
pendent equations. In fact, the existence of perfect MMES for n < 6, n ^ 4, 
will be proven by explicit construction in Section [5] while it is known [45l 
42,43,44] that for n > 8 they do not exist. The case n = 7 remains open, 
although there is numerical evidence that no perfect MMES exist [3T]. In 
conclusion, apart from some special small values of n, not all bipartitions 
can have minimal purity (maximal entanglement) and the requirement that 
a given bipartition be in a maximally entangled state collides with the same 
requirement for a different bipartition. Thus, the bipartitions of a general 
MMES are in a frustrated configuration, and this makes the whole subject 
richer and very interesting. 

Since, according to the structure theorem [15j a perfect MMES can have 
a uniform population probability vector (|zfc| 2 ) = (1/N,--- ,1/N) with 
N = 2™, in Section [5] we focus on this class of uniform states, and restrict 
our quest to it. We will explicitly construct perfect MMES with uniform 
population, and will easily characterize the fully factorized states, i.e. the 
maximizers of the potential of multipartite entanglement that have uniform 
probability vectors. By pushing even further our simplifying assumptions, 
we will explicitly show that, at least for n < 6, the potential admits min- 
imizers and maximizers in the very restricted class of uniform states with 
real phases, (zk) — (Cfc/\/^V), with £fc e {+1,-1}. This states can be nat- 
urally mapped onto the set of binary sequences of length N = 2™, and the 
potential of multipartite entanglement becomes a quartic Hamiltonian on 
binary sequences (or classical spins). It is then quite natural to investigate 
whether there is any relation between the minimizing sequences of 7tme and 
the low correlation sequences that minimize similar long-range Hamiltonians 
studied in [34,35,9,10], which quantify all possible correlations in a binary 
string. However, we will leave this problem for a future publication. 
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A final remark is in order. The study of the minimizers of 7r ME (z) can 
be embedded in a statistical mechanical framework [20] , Let us consider 
the partition function of a system with Hamiltonian ttme(z) at a fictitious 
temperature P^ 1 , 



where /i is the uniform measure (of typical states) on the hypersphere 
{zeC N \ J2 k \ z k\ 2 = 1} induced by the Haar measure on U (N) . The value 
of the free energy Fn(/3) = — lnZ/v, will become that of the minimum 
of the Hamiltonian ttme(z) when the temperature tends to zero, that is 
(3 — > oo, and only those configurations that minimize the Hamiltonian sur- 
vive, namely the maximally multipartite entangled states. In general /3, as 
a Lagrange multiplier, fixes the average value of entanglement, larger values 
of (3 corresponding to a higher multipartite entanglement. In particular, for 
[3 — > one is looking at the typical states. Remarkably, there is a physi- 
cally appealing interpretation even for negative temperatures: for f3 — * — oo, 
those configurations are selected that maximize the Hamiltonian, that is 
fully factorized states. 

This approach has proven to be very useful in the (much simpler) case 
of bipartite entanglement, when the potential of multipartite entanglement 
reduces to the purity ~k a of one of the two component subsystems, and in 
the thermodynamical limit N — > oo, the existence of two phase transitions, 
characterized by different spectra of the reduced density matrices, has been 
shown [23] . 

In order to investigate the statistical mechanics of the richer and more 
complex case of multipartite entanglement, and possibly to unveil its phase 
transitions, it is necessary to study in detail the structure of the potential 
of multipartite entanglement 7Tme and of the highly entangled states that 
give rise to its low energy landscape, the MMES. This paper is completely 
devoted to such a study. 

2 Bipartite entanglement 

In this section we will set up the notation and we will prove some results 
about bipartite entanglement that will be used in the following. We will 
show how the entanglement of a bipartite system in a pure state is related 
to the non- vanishing eigenvalues of the reduced density matrix of one of its 
parts. In particular, in Corollary [3] we will show that in an unentangled, 
separable, state of a composed system the reduced density matrices of its 
two parts are pure, i.e. are 1-dimensional projections and thus have only 
one non- vanishing eigenvalue, that equals 1. On the other hand, a bipartite 
system is in a maximally entangled state if and only if the reduced density 
matrix of its smaller part is completely mixed, i.e. is proportional to the 
identity operator and all its eigenvalues are equal and different from 0. 




6 



Paolo Facchi 



Therefore, as a measure of bipartite entanglement one can use the pu- 
rity, i.e. the sum of the squared eigenvalues, of the reduced density matrix 
of the smaller party. We will do this in Definition [8j One can show that 
purity ranges in a compact interval, its minimum corresponding to maxi- 
mally bipartite entangled states and its maximum to the bipartite separable 
ones. This simple result, which is the content of Lemma [3j together with 
the explicit expression of the purity as a function of the Fourier coefficients 
of the state, as given in Theorem [6] and its corollary, will play a crucial role 
in the following. 

Let us start with some basic definitions. 

Definition 1 (Qubit) A qubit (or spin,) is a quantum system with a two- 
dimensional Hilbert space f) = C 2 . The computational basis {|0),|1)} is a 
privileged orthonormal basis. 

Definition 2 (System of qubits) A system S = {1, 2, . . . , n} of n qubits 
is a quantum system with a 2™ dimensional Hilbert space Hs = (^igsbi; 
with t)i = C 2 . Its pure states are the normalized vectors \ip) € Hs with 
(ip\ip) = 1, and can be expressed in the computational bases as 

\ib) = J2 z k\k), z k eC, M 2 = 1 > (!) 

where k = (fci)ies = (ki, • ■ ■ , k n ), with ki G X = {0, 1}, and 

\k) = \k) s = (g)\k) u {ki^eh. (2) 

ies 

Definition 3 (Bipartition) A bipartition of the system S is a pair (A, A), 
with 1 < ha < n^, where A C S, A = S\A (i.e. S = A + A) and ha = \A\. 
the cardinality of A. The bipartition is said to be balanced if A is maximal, 
that is ha = [n/2] (and ha = [{n + l)/2\), with [x] — integer part of x. 

Remark 1 There is a one to one correspondence among bipartitions and 
nonempty subsets of S of dimension not exceeding n/2. Given a bipartition 
(A, A), the total Hilbert space is accordingly partitioned into Tts = TLa ® 
Ha, where Ha = ®i£A §h w ^h Na = dim 7^ = 2" A , is the Hilbert space 
of the ensemble A of ua qubits. 

Definition 4 (Entanglement) A state \tp) £ Hs is said to be separable 
with respect to the bipartition (A, A) if it can be expressed as a tensor product 
= \4>)a <8 |x) A f or some \4>)a G Ha and \x)a G T^-A- A state that is not 
separable is called entangled. 

The following lemma is a powerful tool in the study of entanglement. 

Lemma 1 (Schmidt decomposition) Given a bipartition (A, A), every 
state \ip) € Hs can be written in the form 



(3) 
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with A fe > 0, £A fc = 1, Y C X nA , and where {\u k ) A } C {|«*>a} c W a 
are orthonormal sets. The set of Schmidt coefficients {A^} is unique. 

Proof In the computational basis — J2iex™a J2 m ex n A Um\l)A ® 
The matrix i = (tj TO ) considered as an operator from to C Wa admits a 
singular value decomposition t = YlkeY V~^k u v 7 f° r some y C I™ 1 , 
with {w^} C and {«<*)} C orthonormal sets, and A fe > J3DJ. 
One gets 1 = (?/>|?/>) = Tr (t*t) — J2^k- The desired result immediately 
follows, with \u k ) A = J2i u< i k) \ l )A and \v k ) A = E m c ™ l^i- n 

It follows immediately that 

Theorem 1 A bipartite state \ip) is separable with respect to the bipartition 
(A, A) iff the set of Schmidt coefficients reduces to {1}. □ 

Remark 2 In general, one wants also a measure that quantifies the entan- 
glement of a bipartite state, i.e. how much the given state differs from a 
separable one. To this purpose, note that one can associate to the Schmidt 
coefficients of a given bipartition {A;, I E Y C X nA } a probability distribu- 
tion p over the finite space X nA , in the following way: 

forZer 

(4) 

otherwise. 

Therefore, it is natural to consider as a measure of bipartite entanglement 
the distance of the probability vector p from the set SEP of the separable 
vectors, concentrated at a point, 

SEP = {p(-) = Se(-)}i £X "A, (5) 

where 6e(-) = 61 n(-). Here 6c is the characteristic function of set C, 

6c(x) = { 1 '* X \° (6) 
W [0 if x i C. K 1 

We consider the distance derived from the L 1 norm, 

d(pi,P2) = l £ \pi(l)- P 2(l)\. (7) 
ieX"A 

It is easy to see that < d(p\,p2) < 1 and that d(pi,p2) = J2i[pi(l) — 
p 2 (l)} + , where [•]+ denotes the positive part. Therefore, if p is the probability 
vector associated to the Schmidt coefficients {A;} of the state in the 
bipartition (A, A), one gets 

min d(p,q) = 1 — max A;. (8) 
geSEP 2 




This motivates the following 
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Definition 5 (Entanglement measure) A measure of the entanglement 
of state \ip) with respect to the bipartition (A, ^4) is given by 

^(l^)) = ^ I (l-maxA i V (9) 

where the maximum is taken over the set of the Schmidt coefficients of 
the state in the given bipartition and Na = 2™ A . 

By noting that, due to normalization, iV7 < maxfc A& < 1, it follows that 

Theorem 2 One gets < S A (\ip)) < 1. Moreover £ A {\^)) = iff \ip) is 
separable with respect to the bipartition (A, A) . □ 

On the other hand, states that maximize the entanglement measure £a are 
the main interest of this paper 

Definition 6 (Maximally bipartite entangled states) A state \tp) that 
satisfies £a (\?P)) — 1 is called a maximally bipartite entangled state with 
respect to the bipartition (A, A) . 

Theorem 3 (Local unitary invariance) A state \ip) is maximally entan- 
gled with respect to the bipartition (A, A) iff 

\i>)=N A 1/2 U A \l) A ®U A \l) A , (10) 

where U A and U A are (local) unitary operators in Ha and H-a, respectively. 

Proof A state is maximally entangled iff £a{\iP)) = 1, i.e., maxj. = 
1/Na- Thus its probability vector is completely mixed, A k = I/-/V4 Vfc G 
X nA . From Lemma [T] one gets the thesis where U A (U A ) is the local uni- 
tary operator in TLa 0~La) that transforms the computational basis into the 
Schmidt one, namely U A \l) A = \ui)a {U a \1)a = \ v i)a)- D 

Remark 3 Note that Eq. (fTOjl implicitly assumes an arbitrary embedding of 
X nA in X A and thus, when > ua, relies on an arbitrariness in the choice 
of the subset {|Z) aIigx^a of the computational basis of party Ha- 

Remark 4 Note that, while separable states (0 are associated to extremal 
probability vectors, concentrated at a point, maximally entangled bipartite 
states are associated to completely mixed probability vectors, uniform on 
X nA . By Theorem [21 the above property can be used as an equivalent defi- 
nition of maximally entangled bipartite states. It has the advantage of being 
independent of the particular measure £a- 

An immediate consequence of Theorem [3] is the following 
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Corollary 1 A maximally bipartite entangled state has the following Fourier 
coefficients in the computational basis 

where Na — 2 nA and l/A, = (Ia\U a \1' a ) with U A the local unitary operator 
in TLa that transforms the computational basis into a Schmidt one. □ 

In fact, there is a link between the set of Schmidt coefficients and the reduced 
density matrices of subsystems A and A. Recall that 

Definition 7 (Reduced density matrix) If p is a density matrix on 
TLa ® TLa, then the reduced density matrix pa is a density matrix on TLa 
defined by 

p A = Tr A p, (12) 
where TrA * s the partial trace over subsystem A. 

Remark 5 The reduced density matrix represents the state of a subsystem 
A, since it determines the statistics of every (local) observables of A. 

Lemma 2 (Reduced density matrix eigenvalues) f5]/ Given a pure 
state in TLs, the reduced density matrices pa and pa of subsystems A and A 
have the same eigenvalues and multiplicities, except possibly for the eigen- 
value 0. 

Proof From LemmaU one gets pa — Tr^ \ip) (ip\ = ^ \k\uk)(uk\, and pa = 

Remark 6 From the proof of Lemma[2]one sees that the Schmidt coefficients 
of a bipartite state are the nonzero eigenvalues of the partial density ma- 
trices of the two parties (and the vectors of the Schmidt decomposition are 
the corresponding eigenvectors). Therefore, from Remark[2]and Defmition[5] 
we obtain 

Corollary 2 Given a state \tp) € TLs an d a bipartition (A, A) one gets 
£a(M) = ^^min{||^ - I^IHj : \<f>) G TLa, {M = 1} 

= ^(i-IM), (13) 

where || • ||x = Tr \ ■ \ is the trace norm and \\ ■ \\ is the operator norm. □ 
Moreover, 

Corollary 3 Given a bipartition (A, A), a state \ip) £ TLs is separable iff 
Pa = for some normalized \(j)) £ TLa o,nd is maximally entangled iff 

p A = l/N A . □ 
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As an alternative measure of the bipartite entanglement between the 
two subsets, which is more suitable to analytical treatment, we consider the 
linear entropy of subsystem A. 

Definition 8 (Linear entropy and purity) A measure of the entangle- 
ment of state \tp) with respect to the bipartition (A, A) is given by 

c A (m = 1 ^ Y (i-M\m, (14) 

where Na — 2 nA , and 

tta(|V» = TtaPa, Pa = Tr A (15) 
is the purity of subsystem A. 

By noting that ita{\4>)) — Y^itfi where {A/} is the set of the Schmidt 
coefficients of the state in the given bipartition it follows that 

Lemma 3 (Purity bounds) Given a state \ip) S Hs and a bipartition 
(A, A), one has 7Ta(|^)) = t^a(\'4'}) an d 

l/N A < Tr A {W) < 1. (16) 

Moreover, 7ta(|t/>)) = 1 and tt aUtP}) = 1/Na iff\4>) is, respectively, separable 
and maximally entangled with respect to the given bipartition. 

Proof The quadratic form ita{\4>)) = Xf reaches its extremal values in 
the simplex A Na — {(A;)/ s x"a|0 < A/ < l,Yli^l = !}• The maximum is 
reached on the frontier, A; = 5g(l) for some 1 < I < ha, while the minimum 
is attained at the interior point where dwA(\ip)) — 0, i.e. A; = 1/Na- By 
Theorems [T] and [3] one gets the thesis. □ 

It follows that Ca^iP)) has a behavior similar to £^(|-0)). In particular, 

Theorem 4 (Linear entropy bounds) One gets < Ca^)) < 1- More- 
over Ca(\iP)) = iff \ip) is separable with respect to the bipartition (A, A), 
while CaH^P)) = 1 iff \ip) is maximally entangled with respect to the bipar- 
tition (A, A). □ 

Remark 7 Let us consider a system composed of an even number n of qubits 
and a balanced bipartition (A, A). The information contained in a maxi- 
mally bipartite entangled state \ip) is not locally accessible by party A or A, 
because, by Corollary[3l their partial density matrices are maximally mixed, 
Pa — Pa — V-^A- Rather, all information is totally shared by them. Note 
that if n is odd, according to Lemma [21 pa cannot be maximally mixed. 
Rather, pa — P/Na, where P = 1 — \v){v\ is a codimension-1 projection, \v) 
being the normalized eigenvector belonging to the eigenvalue 0. Note that 
it is the constraint that the total system is in a pure state that prevents pa 
from being of maximal rank. 
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If the bipartition is not balanced, one gets 

Theorem 5 (Smaller subsystems) A state maximally entangled with 
respect to the bipartition {A, A), is maximally entangled with respect to every 
bipartitions (B, B) with B C A. 

Proof The Theorem is a consequence of Corollary [3] and the property that 
if subsystem A has a maximally mixed density matrix, p A — 1/N A , the 
density matrix of every smaller part B C A is again maximally mixed, 
Pb = Tr^nA Pa = 1/N B . □ 

The explicit expressions of the reduced density matrix and its purity in 
terms of the Fourier coefficient of the state are given by the following 

Theorem 6 (Fourier expression of purity. Form 1) Given a bipartition 
(A, A) and a state \ip) G Hs, one gets 

PA= ^2 z kZlSk A .l A \k A )(lA\ (17) 
k,l£X n 

and 

ka{z) = z k z k 'Ziz v 5 kArA 5 k ^ lA 5k A ,iJk' A ,i' A , (18) 

k,k',l,l'eX n 

where k A = (h) ieA , \l) A = <& ieA \h)i £ H A , 5 k ,i = Si <k = 6{ k }{l), and 
^Ww^S 2 " S 2 ^ 1 = {z e : 5> & | 2 = 1}, (19) 

with N = 2™, are the Fourier coefficients of \if) in the computational basis, 
introduced in Definition^ 

Proof State \ip) can be written accordingly to the bipartition (A, A) as 

W) = Z k\ k A)A ® \k A ) A , 

k£X" 

By plugging this expression into that of p A and ir A given in Definition [8] the 
results follow. □ 

Remark 8 Consider a reference bipartition into two blocks of contiguous 
qubits (C, C), namely C = {1, 2, . . . , n A }, then 

7rc '( Z ) = ^ ^ z (Z,m)^(i',m')%',m)%m')i ( 2 0) 

l,l'eX n A m ,m'eX n A 

where (I, m) = (h, . . . , l nA , mi, . . . ,m„ A ) G X". 

Note that A = p(C) for a suitable permutation p of 5. In fact, there is 
a bijection, 

$: P eCK(p(C), P (C)), (21) 
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between the subset 

V^ A ={p&Vn\p{i)<p(i + l),l<i<n-l,i^n A } (22) 

of the permutation group V n and the set of all bipartitions (A, A) of dimen- 
sion riA- We can write 



(C))p{C) 

= Yl E z p- 1 ((i,m))IOc® \m)c, (23) 
whence, for A = p(C), 

Ka(z) = ^ ^ 2 ; p- 1 (i,-m) ;2: p- 1 (i',m')^- 1 (i',m)^p- 1 (i,m')- ( 24 ) 

For generic bipartitions we have the following 
Corollary 4 (Purity. Form 2) 

= E E E Zfc Zfc © i © m ^ fe ©' ^ESm (^S) 

fcGX s ;ex A me x^ 

were X A and X A are viewed as subspaces of X s with the natural injection, 
and a b — (at © &i)igs = (a; + bi mod 2)i £ s is i/ie XOR operation. 

Proof By substituting in (fTS)) k' = h ® k, one gets 

7rA ( z ) = E z kZk(BhZizi'Sk A ,i' A Sk A <sh A ,i A Sk A ,i A Sk A S)h A ,i' A 

k,h,l,l'eX n 

= E z kZk@hZlZl'0~k®h A ,l&k®h A ,l' 



k,h,l,l'eX" 

z kZk(BhZk®h A Zk®h A1 

k,heX" 



which is the first desired equality. The second equality follows by the iden- 

A 



tifications I = Ha G X" 4 and to = hx € -X'' 4 - □ 



Remark 9 The space X™ is an n-dimensional vector space over the finite 
field X = 7Li with the standard addition and multiplication mod 2. In this 
respect the XOR operation is nothing but the usual sum of vectors of X s 
and X A and X A are vector subspaces. 
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Remark 10 Note that (|25|) can be split into three parts 

+ Yj ^ t Zfe Zfc ffi'® m ^fcem] ' ( 2 ^) 

i£X s iGX^ m€X} 

where = X A \{0}. 

It is an easy exercise to check that the number of monomials |zfc| 4 is 
AQot = 2", the number of monomials |zfc| 2 |z/i| 2 with k ^= h is 

A/tot = 2" (2™- 4 + 2 n * - 2) , (27) 

and the number of monomials Ke[zkZiz m z n ] with distinct indices is 

A/^t = 2™ (2 nA - 1) (2™-* - 1) . (28) 

One gets A/^t + -M-ot +-^tot = 2 2n , in agreement with the first equality in 
(|25p . Moreover, the number of distinct monomials of the various types are 
AA« =A/£ t \ AA( 2 ) = A/^A and ^ =A^/4. 



3 Multipartite entanglement 

The aim of this section is to generalize the ideas of the previous section to 
the case of multipartite entanglement. We require that the information in 
a maximally multipartite entangled state be distributed as well as possible. 
In the ideal case this would mean that 

Definition 9 (Perfect MMES) A state maximally entangled with re- 
spect to every bipartition (A, A) is called a perfect maximally multipartite 
entangled state (perfect MMES) . 

Theorem 7 (Perfect MMES characterization) The following state- 
ments are equivalent: 

1. \ip) £Hs is a perfect MMES; 

2. pa — for every subsystem A C S with ua < n/2; 

3. pa — I/Na for every maximal subsystem A C S; 

4- £a(\i/>}) = £a{\iI>)) — 1 for every balanced bipartition (A, A); 
5. tta(\iP)) — for every balanced bipartition (A, A). 

Proof Equivalence between 1 and 2 follows from Definition [5] and Corollary 
[21 Statements 2 and 3 are equivalent by Theorem [5j Equivalence between 
3 and 4 follows from Definition [6] and Theorem |4j Finally, 4 and 5 are 
equivalent by virtue of Lemma G3 □ 
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Remark 11 Note that the requirement that a given balanced bipartition 
(A, A) be in a maximally entangled state could collide with the same re- 
quirement for a different balanced bipartition (B,B), with B ^ A. Indeed, 
the local unitaries U A and U A in Theorem [3] are in general nonlocal for 
the bipartition (B,B). Thus, at variance with the bipartite perfect 
MMES cannot exist. This motivates the following definition. 

Definition 10 (MMES) Let us define the potential of multipartite entan- 
glement as 

W|^>)=(r/ 21 ) E MM)' (29) 

|A| = [n/2] 

A maximally multipartite entangled state (MMES) \(p) is a minimizer of 

7TME, 

7r ME (|¥>)) = irrm{7T ME (|V>)) I W G H 8 , = 1}- (30) 

The potential 7Tme measures the average bipartite entanglement over all 
possible ( r ) balanced bipartition and thus inherits the bounds on the 
purity given in Lemma [3j namely, 

Lemma 4 (Bounds on 7Tme) The potential of multipartite entanglement 
satisfies 

l/N A < 7r ME (|V)) < 1, (31) 

with N A = 2l"/ 2 l ; for all normalized \ip) G Hs- □ 

The upper and lower bounds are characterized by the following 

Theorem 8 (Optimizing states) The upper bound ttme(z) — 1 is at- 
tained by the fully factorized states, whose Fourier coefficients in the com- 
putational basis z = (zk)keX n are z k — YiieS a k > ™^ l a ol 2 + \ a \\ 2 = 1- 
On the other hand, the lower bound ttme(^) = */ attained, would 

correspond to a perfect MMES. 

Proof 7Tme(|V')) — 1 iff 7r A(|V ; )) = 1 f° r an balanced bipartitions (A, A). 
By Lemma [3] this happens iff is separable with respect to all balanced 
bipartitions. Now, note that 1^) = \vi)a ® \v2) a an d \ip) = \vs)b ® \ v i)s 
iff \ip) = |ui3)Ans ® |«i4>Anfl ® \v2z) Ar\B ® \ v 2i) AnS- Since for au * G S 1 , 
{z} = n r A r for a suitable set {A r } of maximal subsystems, one has \ip) — 

0ies \ v i)i witn = L Thus = S feeX " FliGS^^I^)' and the first 

part of the theorem follows by setting a\. = Concerning the second 

part, 7Tme(|V')) — iff tta(IV')) — for all balanced bipartitions 

(A, A). By Theorem [7] this happens iff \ip) is a perfect MMES. □ 

Remark 12 In words, a perfect MMES is characterized by a multipartite 
entanglement that is maximum, in the sense that it saturates the minimum 
of the purity and such a minimum does not depend on the bipartition. 
However, if the minimum of the potential of multipartite entanglement is 
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strictly larger than the lower bound in Lemma 21 i.e. min7TME > 1/N A , it 
may happen that different bipartitions yield different values of n A , some of 
them smaller than min7TME, some larger. In such a situation, one can seek 
those states among the minimizers that have the lowest variance. This quest 
can be recast as an optimization problem [21] . We will not elaborate further 
on this issue. 

Now we will examine in more details the potential of multipartite entangle- 
ment and we will determine its form. 

Theorem 9 (Fourier expression of 7Tme- Form 1) Given a state \ip) 6 
H.S, the potential of multipartite entanglement has the following expression 
in terms of its Fourier coefficients in the computational basis z = (zfc)fcex™ 

kme(z) = A(k,k';l,l';[n/2])z k Zk' ziz v , (32) 

with a coupling function 

A(k, k'; I, I'; n A ) = ^A(k, fc'; I, I'; n A ) + ^A(k', fc; I, I'; n A ). (33) 

where 

A(k,k';l,l';n A )=(™) £ S^J^Ju^S^. (34) 

V J \A\=n A 

Proof The result follows by plugging the expression (|18|) of ir A given by 
Theorem [6] into (|29|) of Definition I10| and by symmetrizing. □ 

Remark 13 In the spirit of Remark it is easy to see that the average can 
be extended to the whole permutation group, yielding 

A(k,k';l,l';\C\) = — St. 6h> i ,„><5fc i .^.StJ v .(35) 

Remark 14 Note that A would have served as well as A as a coupling func- 
tion, namely 

ttme(^) = ^2 A(k,k';l,l';[n/2]) z k z k ' zizi> . (36) 

k,k',l,l'eX n 

However, while 

A(k, fc'; I, I'; n A ) — A(l, I ; fc, k ; n A ), 
which ensures the reality of 7Tme7 one gets 

A{k , fc; I, I'; n A ) = A(k, fc'; I , I; n A ) — A(k, fc'; /', I; n A ). 

Thus, A{k, fc'; I, I'; [n/2]) is a symmetric function of the pairs (fc, fc') and 
(1,1') only if n is even, when [n/2] = n/2 = n A = n A . Since 7Tme does not 
depend on the antisymmetric part of the coupling function, we shall use the 
symmetric coupling function A. We summarize its properties, which easily 
derive from this Remark in the following 
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Lemma 5 (Coupling function symmetries) The coupling function A : 
X 2n x X 2n xN^Q has the following symmetries 

A(k, k'; I, I'; n A ) = A(k', k; I, I'; n A ) = A(l, I'; k, k'; n A ), (37) 

for every k, k' , I, I' € X n and every n A with 1 < n A < n — 1. □ 

The following definition and the subsequent lemma are the main ingredients 
for determining the explicit expression of the coupling function A. 

Definition 11 (Admissible set) Let us define the admissible set as the 
set of all quadruples of sequences that yield a nonvanishing contribution to 
the function A, that is 

Qn A = {(&! G X in \k A = l' A 7 k' A = l A , and k A = l A , k' A = l' A , 

for some (A, A) with \A\ = n A }. (38) 

Obviously, 

Qn A cQ= |J Q s . (39) 

0<s<n 

Lemma 6 (Admissible set characterization) The set Q is the kernel 
of the function q : X 4n — > X n , 

q (k, k', i, i') = ({k © o v (k' © i')) a ((fc e i') v (k' © 0) , (40) 

where a © b = (dj © &i)igs = (dj + bj mod 2)j £ s is t/ie XOR operation, 
a V = (aj V bi)i e s = (di + bi + aibi mod 2)j £ s i/ie OR operation, and 
a A b = (a,i A bi) ie s — {aibi)ies the AND operation. 

Proof The proof consists in a straightforward application of the above de- 
fined binary operations: 

Q = {(k,k',l,l')\k A = l' A , k' A = I A , k A = l A , k' A = l' A , for some A C S} 
= {h = k[ = li, kj = lj, k'j = with % e A, j e A} 
= {h © I'i = 0, k'i © h = 0, kj © Z,- = 0, k'j © with i e A, j e A} 

= {(fc 4 © /D v (fc( © ij) = o, (kj © y v (fcj © l'j) = o, * e A j e A} 
= {((** © O v (fc- © ii)) a © /,) v © i'i)) = o, i e S} 

= {({k © J') V (fc' © 0) A ((fc © V (fc' © I')) = 0} 
= kerq. 

□ 

Remark 15 Note that X n can be viewed as a product ring (of n copies of X = 
Z 2 ) with the addition and multiplication mod 2 defined componentwise, as 
usual. In this respect, the XOR operation is the sum a + b and the AND 
operation is the product a ■ b of elements a and b of the product ring X n . 
The OR operation is nothing but a + b + a ■ b. 



Multipartite entanglement in qubit systems 



17 



After having proven all preparatory lemmata, now we come to the main 
result of this section that establishes an explicit form for the coupling func- 
tion of the potential of multipartite entanglement. 

Theorem 10 (Coupling function) The coupling function A has the fol- 
lowing expression 

A(k, k'; 1, 1'; n A ) = g((k © Z) V [k' © I'), {k © I') V (k' © l);n A ) , (41) 
where 

g(a,b;n A ) = S (a A b) g(\a\, \b\; n A ), (42) 
with \a\ — ^2 ieS Oi, and 

" 9 ^ nA) = \{n A ) 
Proof Let 

a = (fcffi I) V [k' © V) and b = (k © /') V (Ar' © I). 

By Lemma [6] (k, k', I, I') £ Q iff a A b = 0, and thus Si = {i € S\ai = bi = 
1} = 0. Therefore, 

S = S +A 1 +B 1 , 

where 5 = {i € SK; = 6j = 0}, A\ = {i € £|<Xj = 1}, and Bi = {i £ 516, = 
1}. Moreover, it is easy to see that a,; = iff ki = h and k\ = l[, with i £ S. 
Thus, if (k,k',l,l') £ Q, then fc^ = Z^ and k'^ — l'j, and, analogously, 
kg i — I'g and fc^ = l§ 1 . As a consequence, (k, k' , I, I') £ Q nA iff there is a 
repartition (A, A), with |A| = n A , such that 

icii and A C -Bi , 

that is ^! C 4 and £?x C A In other words, (fc, k', I, V) £ Q„ A iff (k, k', I, V) £ 
Q and |Ai| = \a\ < n A , \B\\ = \b\ < n A . Therefore, we can write 

(k, k', I, I') G Qn A iff a A b — 0, with \a\ < n A , \b\ < n A , 

whence 

A(k, k'; I, I'; n A ) = 5 (a A b) <%,, M (|a|) 5 [0 , nj[] (\b\) ( ^ ) #(*, I, I'), 

where #(fc, k', I, I') is the number of terms of the sum (|34[) that contribute 
to the function A in Theorem [5J 

Now, according to the above conclusions, for a given admissible quadru- 
ple (k,k',l,l') £ Qn A t ne number of terms #(k,k' ,1,1') is given by the 
number of bipartitions (A, A) with |A| = n A and with A C -Bi = A\ + So 
and A C Ai = Bi + So- Since A n A = 0, Ai C A and B\ C A, parties A 
and A contend only for So — So A + So n A, namely 

A = A : + S n A and A = B x + S n A. 
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Thus, their number equals the number of ways that |^4\Ai| objects can be 
chosen from among |5o| objects. But = \A\ — \A\\ = n A — \a\ and 

|5 | = \S\ - - = n - \a\ - \b\. Therefore, 



#(k,k',l,l') 



n-\a\-\b\ 
n A — \a\ 



By putting all together, and by stipulating that the binomial coefficient 
is zero when its arguments are negative, we obtain the stated form of the 
functions A and its symmetric part A. □ 

Remark 16 It is not difficult to see that an alternative form of g is the 
following 



g(s,t;n A ) = 



1 

2 \s,t 



n A \ n A 



n A 
t 



n A 
s 



(44) 



where 



s,t) s\t\(n- s-t)\ 
is the multinomial coefficient. 

By using the explicit form of the coupling function A one can give the 
potential of multipartite entanglement a different form that has the advan- 
tage of being a sum over three indices only. 

Theorem 11 (7Tme- Form 2) The potential of multipartite entanglement 
can be written as 

ttmeM = ^2 g{l, m; [n/2]) Re [z k z k(Sl(Bm z km z mm ] . (45) 

Proof Since k © = k, 

Q = {(k, k', I, l')\k A = l' A , k' A = l A , k A = l A , k' A = l' A , for some A <Z S} 
= {k A = l' A , k' A ®l A = 0, k A = l A , k' A ®l' A = 0} 
= {k A = l' A @k' A @ l A , k' A @l A = 0, k A = l A @k' A @ l' A , k' A @l' A = 0} 
= {k = k' © I © I', k' A © l A = 0, k' A © l' A = 0} 
= {k = k' © I © I', (k' © A (k' © I') = 0}. 

Moreover, since k © k = 0, substituting for k = k! © I © V one gets a = 
(k © I) V (k' © V) = k' © I' and b = (k © I') V {k 1 @l) = k'®l. Therefore, 



A(k, k'; I, I'; n A ) = 5 k ,k'®l®l' g{k' © I', k' © I; n A ), 



whence 



ttme(z)- 9(k' ®l',k' ®l;[n/2])z k , i 

k',i,i'ex n 



z k ' zi z v . 
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By setting I' = I® k' and I = m® k', one obtains 

kme{z) = ^2 d(h w; [n/2]) z k / @ i &m z k > z k >® m z k , 



The thesis follows from the reality of ttme(z)- □ 

In analogy with the bipartite case examined in Remark I lOi the sum in (|45 
can be split into three terms. 

Corollary 5 (7Tme- Form 3) 



ttme(*)= £ M 4 + 2]T 9(\l\A[n/2])\z k \ 2 \z k!Sl \ 2 
hex™ fcex™ lex? 

+ X! 9{l,m][n/2])Ke[z k z kmmm z m iz kmm ], (46) 

hex™ i.mex" 



w/iere X^ 1 = X"\{0}. 

Proof The monomials |zfc| 4 are obtained from (|4"5|) when I = m = 0. In such 
a case g(0,0; [n/2]) = g(0,0; [n/2]) = 1. On the other hand, the monomials 
\z k | 2 \zh | 2 with k ^ h are obtained when either I — or m = 0. In such a case, 
since £ (MO) = 1 for all I G X n ,g(l,0; [n/2]) = ff (Q,Z; [n/2]) = g(\l\,Q; [n/2]). 
□ 

A measure of the complexity of the potential of multipartite entanglement 
is given by the number of its terms. In particular, as we will see in the 
following, the crucial ones are the interfering monomials Ke[z k ziz m z n ]. 

Theorem 12 (Number of terms in 7Tme) Consider ttme{z) . The number 
of distinct monomials \z k \ 4 and the number of distinct monomials \z k \ 2 \zh\ 2 
with k ^ h are 



" ~ + 3+(-i)»U"/ 2 ]y ' 1 j 

respectively. The number of distinct monomials Ke[z k ziz m z n ] with distinct 
indices is 



^(4) =2 „-3 E ( n )( n- S)=2n _ 3 E (" .(48) 
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Proof The total number of terms of the sum in (|45jl is given by 
Mot= E S Q.(9{l,rn;[n/2})) 

k,l,m£X n 

= 2" Yl S o(lAm)6 [Qm (\l\)d [om (\m\) 

l,m£X™ 

= 2 " E E S (lAm)S s (\l\)S t (\m\) 

< s ,t<[n+i] i,meX" 

= E E E mm) 



0<s,t< 



i±il Jex* 



-2" E 



< 5 ,t<[i+i] 



n \ n — s 



n \ in 
M 



Therefore, the total number of monomials |zfc| is 
A4V = E Ms(0,G;[n/2]))=2" 

while the total number of monomials |zfc| 2 |zft| 2 with k ^ /i is 
= 2 £ ^ £ Q , (g(l, 0; [n/2])) = 2» +1 £ 

i<*<m 



= 2" Y 



2", 



n \ n 
[t 



n 

n — t 



l<t<[=Ji] v 7 



2 " £ " 



KKn 



2 2n _ 2 n + 2 n 



2^ n 2 n 



[tj<*<« 

+ 2™ 53 

[§]<*<m 

E 

w<*<m 



n 

n/2] 



3 + (-1)" U n / 2 1 



On the other hand, the total number of monomials Ke[zkZiZ m z n ] 
tinct indices reads 



with dis- 



A/7 



(4) 



E E Ms0,m;[n/2])) = 2" E 



1<8,« 2fl 



t 



The results follow, since by symmetry, the numbers of distinct monomials 
are = , = At™ /4, and AA 4 ) = A/$ /«• □ 
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Table 1 Number of monomials in 7tme(«)- 



^(D ^(2) ^(4) 



2 


4 


4 


1 


3 


8 


24 


12 


4 


16 


80 


84 


5 


32 


400 


680 


6 


64 


1312 


4000 


7 


128 


6272 


28672 


8 


256 


20736 


162624 


n — ► oo 


2 n 


22n-2 


2 n-3gn 



Remark i7For large values of n, by making use of Stirling's approximation 
one gets 

n 



n/2 



V TO 



hence, from (|4"T)l 



The asymptotics of A/' 1 - 4 -' is a little more elaborated. First note that, by 
Stirling, 

z — ' - ' - + » * \ \nn 



i<, f< \»+i] limJH- 

where the function 77 : Zi 2 — > R, defined on the simplex Z\ 2 = {(x,y) G 
[0, l] 2 | x + y = 1}, is the entropy 

H(x,y) = -x\ogx - y logy - (1 - x - y) log(l - x - y). 

Then, for n — > oo, by using the same arguments as in the proof of Laplace 
- De Moivre theorem [29], one can show that 

~ 2™- 3 3™, n ^ oo. (50) 

The numbers of different types of monomials appearing in ttme(z), as well 
as their asymptotic expansions, are given in Table Q] 

We conclude the section by exhibiting another form of multipartite en- 
tanglement. 

Theorem 13 (7Tme- Form 4) The potential of multipartite entanglement 
can be written as 

kme(z) = 1 - - y g(l,m;[n/2}) \z k z kS) i @m - z k@ iz k@m \ 2 . (51) 

k,l,m£X n 
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Proof Let us consider (|45|) . By substituting the identity 

1 2 

Re [Zk Zk($lQm ZkQl Zk^m] = \ Zk z k®l®m ~ Zk®l ^fcffim | 

+ 2 \ z k z k®l®m\ + ^ \ z k®l z k<$m\ , 

one gets 

^me(z) = — - #(Z, m; [n/2]) Zfc©z em - z fcei Zfc© m | 2 

^fe^fee/ffiml + \ z k®l z k@m | ) 

Now, by simple manipulations, 

|2 



and 



\ 2 = ^2 \zkZk®i\ 2 ^2 g(l@m,m;[n/2]). 

k,l,m£X n k,leX" mGX" 

Thus, 

t^me(z) = -- ^2 9{Um;[n/2]) \z k z k ®i® m - z k @i z k @ m \ 2 

k,l,m£X n 

+ X! \ z k z mi\ 2 ^2 9( l ® m , f n;[n/2}). 



I 2 

k,leX n m£X 

Let us assume for a moment that 



J2 g(l®m,m;[n/2]) = 1, \JleX n . (52) 



Then, the result follows by normalization (fTi?)) . z € S 2JV_1 , since 



\ z k z k ®i\ 2 = ( J2 \z k \ 2 ) =1 



In fact, equality (|52|) . is a consequence of the following lemma, for tla = 
[n/2]. □ 

Lemma 7 XTie following equality holds 

Y = H fl(Iem ) m;n A ) = l 1 V/ e X™, Vn A € S. (53) 
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Proof From Eq. (|42| in Theorem \W\ we get 

Y= ^2 S ({l®m) Am) g(\l &> m\,\m\;n A ). 



Let us define the set B = {i g 5 | Zj = 0} c S, so that Z^ = 0. We get 
(I (Bin) Am = iff tub = and (lg(Bm^)Am§ — 0. But the second equality 
is identically satisfied, because l B is a vector of all 1. Thus So((l(Bm)Am) = 
So{mg) and we get 



Y = E 5 o(m s )g(\l s ®m s \ 1 \m s \]n A ) 

meX" 

= E fl(|ll...lem| s |m|;n i4 ) 
= E - H> H;tia) 

= E E MH^O-M;^) = E 



meX' 0<i<r 



0<t<r 



<?(r - t,t;n A ), 



where r = \B\ 
Remark [TBI 



Let us now use the form (144|) of the function 5 given in 



Y 



- E 

0<*<r 

2 I r 



n 

r-t,t 



n A 
r-t 



n A 
t 



r - t 



E 

0<t<r 



n A 
r - t 



n A 
t 



r - t 



n A 



r - 1 



0<t<r 

By recalling Vandermonde's identity |16j . 



E 



m \ in — m 



since n A + n A = n, we get 
-1 

E 



Y = 



0<t<r 



k-j 



n A \ I n — n A 

r - t 



= 1. 



□ 



Remark 18 Recall that, by Theorem [8j the potential of multipartite entan- 
glement attains its upper bound ttme(z) = 1 on fully factorized states. Thus 
in (|5ip the nonegative sum which is subtracted from unity represents the 
amount of entanglement of 1-0), and MMES are those states that maximizes 
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the distances \zk z&©i® m — Zk®i ^fce™| 2 - I n fact, the average over balanced 



bipartition of the linear entropy (jT4j) yields 

N A 

Cme(z) = j^j(1-t: me (z)), (54) 

with Na — 2l 7V / 2 ]. Thus, apart from a normalization factor, the sum in (|5ip 
is nothing but the average linear entropy. Note that the number of terms 
in the sum is J\f^ given in (|48[). since the terms with I — or m = 
identically vanish. 

Remark 19 In the spirit of the above Remark, one can prove Theorem 1131 
by following a different path. First, one can easily write an expression anal- 
ogous to ([5]) for the purity tta(z) of a given bipartition. Incidentally, this 
would give an explicit expression of Ca(z). Then, one considers the average 
over balanced bipartitions and, by noting that the proofs of Theorems ITUl 
and 1111 do not depend on the particular form of the monomials ZkZiZ m z n , 
that can be replaced by \zkZ% — z m z n \ 2 , one obtains the desired result. By 
comparing the two proofs, since the average of 1 is 1, one can easily distillate 
an alternative combinatoric proof of Vandermonde's identity. 

Example 1 Consider n = 2 qubits. One gets 

7Tme(z) = |z o| 4 + \zoi\ 4 + \zio\ 4 + kn| 4 

+2 (M 2 M 2 + M 2 M 2 + kn| 2 M 2 + kn| 2 M 2 ) 

+4Re(z o2oi^ioZii) 
= 1 - 2|z 00 zii - z Q1 z w \ 2 . (55) 

The first equality follows from CorollaryEl while the second equality derives 
from Theorem [T51 Note that the number of terms AA 1 ) = Af^ — 4 and 
J\f( 4 > = 1 is in agreement with the counting of Theorem [T2l and Remark [l~8l 
See Tabled] 

Example 2 For 3 qubits we will give the potential of multipartite entangle- 
ment in the form 4 of Theorem [T3j 

ttme(z) = 1 — 2 (j z p(ooo) z p(oii) — 2 p(ooi) z p(oio) | 2 
pec 3 

+ l z p(ioo) z p(in) ~ z p(ioi) z p(no) | 2 
+ g l z p(ioo)^p(on) — ^p(ioi)^p(oio)| 2 

+ ~^\ z p(ooo) z p(iii) ~~ z p(ooi) z p(no) | 2 ) j (56) 

where the sum is over the 3 cyclic permutations 

C 3 = { S l Kl,2,3) = (2,3,l), i = 0,l,2} (57) 
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of the qubits S = {1,2,3}. Here p(k) denotes the natural action of the 
permutation group on k, 

T n x X n 3 (p, k) » P (k) = {k p(i) ) ieS G X n . (58) 

In agreement with Remark 1181 the number of distinct terms is 
3 x 4 = 12. See Table Q] 

Now we will focus on the problem of the existence of perfect MMES. 
In particular we will try to construct them by using characterization 2 of 
Theorem [7l It is not obvious that a state with 7Tme = 1/Na exists: in order 
to find a solution one must solve for pa — 1/Na,V (A, A), and this set of 
equations might not admit a solution. 



4 Perfect MMES. Probabilistic approach 

We will look more closely at the equations 

Pa = for every subsystem A C S with \A\ < n/2, (59) 

that, according to Theorem [3 characterize a perfect MMES. Although we 
could consider only maximal subsets A c S, with \A\ = [n/2], it will be 
more convenient to consider also smaller subsets A. 

Let us first consider the diagonal elements in the computational basis 
{\e) A }iex"A C H A - By Eq. (17]) of TheoremEJ one gets 

(£\pa\£) = z k zi5k A .i A 8k A ,i5i A .i = 2J z k zi5 k ,i5k A .,i 

= £ \z k \ 2 S kA/ . (60) 

Therefore, from (|59|) we obtain 

(£\pa\£)= W\ 2 S kA ,t = l/N A , (61) 
keX" 

with N A = 2l A l, W G VA C S, with |A| < n/2. 

Now note that, due to normalization, ^ |^fc| 2 = T we can look at 
(|zfc| 2 )fceX" as a probability vector on the finite space X n of n classical 
bits. In view of this interpretation, we will introduce the 

Definition 12 (Population probability) Given a state G Tis and 
its Fourier coefficients (z k ) ke x n in the computational basis, we define the 
population probability vector in the computational basis, 



Ps(k) = \z k \ 2 , 



(62) 
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as the probability of the binary sequence k = (fcj)igs G X n . Moreover, let 
E[-\ denote the expectation value with respect to Ps, 

E[f(k)}= £ f(k)Ps(k)= J2 Z( fc )M 2 > ( 63 ) 
hex™ hex™ 

for any function f : X n — > C. 

According to the above definition, Eq. (|61[) reads 

E[S kA ,t] = 2~\ A \ WeI |A| , VA c S, with |A| < n/2. (64) 

By noting that the above expectation value is nothing but the marginal 
probability distribution 

E[S kAA ]=P A (e), (65) 
we have arrived at the following 

Theorem 14 A necessary condition for a state \ip) £ TLs to be a perfect 
MMES is that all the marginals over ha < n/2 variables of its population 
probability vector in the computational basis Ps(k) — \z k \ 2 , are completely 
random: 

P A {l) = 2^ A \ V£e J |A| , VA c S, with \A\ < n/2. (66) 

□ 

Remark 20 For A — p(C) with p £?„ and G = {1, 2, . . . , ua} one can write 



E 



W e A" A , Vp € P n , with n A < n/2, (67) 



which means that 

P p{C ) (£) = 2-' nA , W e , e P„, with n A < n/2. (68) 

Remark 21 According to Theorem 1 141 a first step in the problem of seeking 
perfect MMES is the following: Search for all probability functions on X n , 
whose marginals on ua < n/2 variables are uniform. 

The solution to this problem is given by the following 

Theorem 15 (Perfect MMES population) The population probability 
vector in the computational basis of a perfect MMES of n qubits has the 
form 

\z k \ 2 = p s (k) = 2- n + £ 4 3 II ( 2fc i« - !). keS > ( 69 ) 

%<r<nj£[S r ] l<l<r 

for some cy € R, where [S r ] = {(ji, . • • , j r ) G S^Ui < j 2 < • ■ ■ < jr} 
denotes the set of ordered vectors of S r . 
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Proof Note that any function on X n is a multilinear function of the com- 
ponents of k G X n , because kf = k{. Therefore, we can write 

^(fc)= c (o) +£ e 4 r) n t 2 ^- 1 )' 

r£Sj£[S r ] l<l<r 
(r) 

which depends on the real parameters Cj e R, whose number is 
The normalization of Ps implies that 

i = e PsW = 2" C (°) + E E 4 r) II E ^ - 1) = 2 " c(0) ' 



that is 



c (0) = 2 -n_ 



Let us now consider a subset with one element A — {j}, with j e S. For 
any kj € X one must have 

\ = P{j}(kj) = E P ^(fc) - 2- 1 C (°) + cf\2kj ~l) = \ + cf(2k 3 - 1), 
that is 



k£X* 



4 1} =0, jeS. 



Analogously, for a subset with two elements A — {ji,j2} 5 

2- 2 = P A (k h ,k ja )= E ^(fc) = 2- 2 + cg i ) j2) (2fc J1 -l)(2fc, 2 -l), 



that is 



By induction we get 



4 2) =o, je[n 



c^ r) =0, y.je[S r ], forl<r<n/2, 



and the result follows. □ 



(r) 

Remark 22 The range of the free parameters c!j is determined by the in- 
equalities < Ps{k) < 1, Vfc e S. Their number is given by 



E i[sii= E 

^ <r<n 

4e 



0<r<4 



1 E 

Y < r < n 



= 2 n-l _ ! + (-!)" 



r <r<n 

n 
[n/2] 



n 

n — r 



(70) 
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(r) 

The particular solution cj- = for all r G 5 that yields a uniform proba- 
bility Ps(k) = 2~™ will play a role in the following. 

Theorem [15] completely determines the structure of the moduli — 
\ z k\ = \J Ps(k) of the Fourier coefficients Zk of a perfect MMES in the 
computational basis. However this is only half of the work. In fact, the easy 
one. It remains to determine the phases, defined in the following 

Definition 13 A state G Tls can be expressed in the computational basis 
as 

= X] Zfe l fc )' z k = r k (k, (71) 

where the Fourier moduli belongs to the intersection of the positive hyper- 
octant with a hypersphere 

re (R+) N n^- 1 = {(r k ) keX ~\r k £R+,Y, r t = 1}> N = 2 "> (72) 

k 

while the Fourier phases belongs to the torus T N = (S 1 )^ 

( g T N = {(a) fce x- 1 Cfc e C, |Cfe| = i}. (73) 

We will now show that the phases £ of a perfect MMES are solutions to the 
system of the off-diagonal elements of the equation pa = 1/Na- 

Theorem 16 (Perfect MMES phases) A state G Us is a perfect 
MMES iff its Fourier phases £ in the computational basis are solutions to 
the equations 

m£X A 

V£,£' G X A , i ^ e, VA c S, \A\ = [n/2], (74) 



where r k = W Ps(k), with the population probability vector Ps(k) given in 
Theorem \15l for some coefficients c*p G M. 

Proof The off-diagonal elements of pa in Eq. (fT"7|) of Theorem [6] read 

{Z\Pa\Z') = X ZkZiSk A ,i A Sk A ,i5i At if = X zt®kx^t'®k A ) 

k,leX" k A £X A 

Ve J,' G X A , I ^ £'. Thus, by Equation (jHSJ, Definition [TJ and Theorem [T5] 
one gets the desired result. □ 

Remark 23 An alternative form of (|74[) in terms of the permutation group 
is 

e vp S (k)p S (i)CkCi n Ku^w n n w;=°> 

fc,ZeX" Tt<j<n !<i<f !<J<f 

VpGP n , Wi'el ni , l<n A <n/2. 

(75) 
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Table 2 Number of equations vs number of variables 



n m e m x 



2 4 5 

3 6 12 

4 72 21 

5 120 48 

6 1120 86 

7 1960 192 

8 16800 349 



n^oo ^2/mr 2 2 " 3 2' 



Let us now investigate whether the system of equations (|74[) admits a so- 
lution or not. In particular, it is important to count the number of equa- 
tions and of variables and to look for which values of n the system is over- 
determined. 

Theorem 17 (Number of equations and variables) The set of equa- 
tions \7Jj\j determining a perfect MMES is a system of 



me = 2 [ n /2] (2 [n/2]_ 1) ^ 

real equations involving 



n 

n/2] 



(76) 



m x = 3 • 2«- ^ ( ) (77) 



,n-i l + ("l) n ( n 



real variables. 



Proof By noting that exchanging I and I' one obtains the complex con- 
jugate, the counting of real equations coincides with the total counting of 
equations (|74| . Since I ^ we get 

m e = \X A \ (\X A \ - 1) #(A), 
where #(A) is the number of maximal subsets A <Z S. Now, \X A \ — 2l"/ 2 l 



and #(^4) = I j /2] y ' an< ^ follows. On the other hand, the variables 

are the 2™ phases £ and the parameters , whose number is given by (fT0|) , 
for a total number of m x variables. 

Remark 24 For large values of n, by Stirling's approximation one gets 

m e ~2 2 "W— , m a; ~3 2™- 1 , n -► oo. (78) 
V mr 

As shown in Tabled for n > 4 the number of equations is larger than the 
number of variables and the system is overdetermined. Therefore, symme- 
tries must play a crucial role in order to assure the existence of a solution. 
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4-.1 Examples 

4-1.1 Two qubits Let us consider the case of n — 2 qubits. S = {1, 2} and 
we get from (p9")) 

r 2 k = P {1 , 2} (k) = j(l+ca 1 a 2 ), (79) 

with Ui = (2ki — 1), i € S. Normalization and positivity, < Pi t t 2\(k) < 1, 
Vfc € X 2 imply that c 6 [—1,1]. Equation ([74)) particularizes to 

I JWioCooCio + nnXiiCoiCu = o 
{^oonnCooCoi + riomCioCn = 0, 

and, by noting that Tq = r\ x = (1 + c)/4 and j"q X = r 2 = (1 — c)/4, one 
gets 



J Vl~? (CooCio + CoiCii) =o 

\VT^ ^ (CooCoi + CioCn) =o 

The above system reduces to a single equation 



(81) 



Vl-c 2 (CoodiCioCoi + 1) - 0. (82) 

1. A first class of solutions is \c\ = 1 and arbitrary phases. This yields, for 
c = 1, r i = r 10 = and r Q0 = r n = l/\/2, whence 

|V)=-^(Coo|00)+Cn|ll>), (83) 

while, for c = — 1, roo = ?'n = and roi = ?*io = l/v2 whence 

|^ = ^=(C 01 |01)+Cio|10)). (84) 

The above states are known as Bell states. They are, obviously, maxi- 
mally bipartite entangled. Indeed, for n — 2 multipartite entanglement 
reduces to bipartite entanglement. 

2. On the other hand, when |c| ^ 1, the perfect MMES are 

= ^v / TT^(Coo|oo> + Cn|ii)) + ^\/T^(Coi|oi> + Cio|io)), (85) 

where the phases must satisfy the condition 

Coodi = -CoiCio- (86) 

Therefore, 



ty) = (/J|00> + £|11» - icx^P- ( 7 |01) + 7|10» , (87) 



with a = (CooCii) 172 , = (CooCn) 1/2 and 7 = (CoiCio) 17 



2 
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3. The particular case c = corresponds to a uniform amplitude distri- 
bution rfc = 1/2, k € X 2 . To such a class belong perfect MMES with 
phases that are ±1. 



\ E Ck|fc>, C*e{-i,+i}, n Ob = -i. 



(88) 



kex 2 



hex 2 



4-1.2 Three qubits Let us consider the case of n — 3 qubits. S = {1,2,3} 
and we get from ([69]) 

t\ = Pl,2,a(k) = - (1 + CiCr 2 (T3 + C2Crier 3 + C3Cri(72 + rfcricr 2 cr 3 ) , (89) 

where Oi = (2fcj — 1), with i £ S. On the other hand, from (|74[) we obtain 

(90) 



^000^100 + ^001^101 + ^010^110 + ^011^111 = 
■SOOO^OIO + ^001^011 + ^100^110 + ^101^111 = 
z 000^001 + ^010^011 + ^lOO^lOl + ^110^111 = 



Note that the three equations are obtained by a cyclic permutation of the 
three qubits S. 



1. If Ci = (i € 5) and d £ [—1,1], one gets from 



1 m j \ i-(-i) |fc| MI 

- (1 + da 1( 7 2( 73) - / ' . 



lei 3 



Thus, 



Vl — d 2 (CoooCioo + C001C101 + CoioCno + CoiiCm) = 
Vl — d 2 (CoooCoio + CooiCon + CiooCno + CioiCm) = 
^ Vl — d 2 (CoooCooi + CoioCon + CiooCioi + CiioCm) = 



(a) If \d\ = 1 the phases are arbitrary and the MMES is 



\1>) = 2 (Cooi|001> + Coio|010) +Cioo|100> + Cin|Hl>) 



= 5 E Cfclfc), 



|fe| odd 

when d = 1, and 

W =^(Cooo|000) + Cqh|011)+Cioi|101) + Ciio|H0)) 



= \ E c*i*>, 



(91) 



(92) 



(93) 



(94) 



|fe|< 



when d 
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(b) When \d\ ^ 1, the phases must satisfy 






C000C100 + CooiCioi 


= a 




C000C010 + CooiCoii 


= p 


< 


C000C001 + CoioCoii 


= 7 




CoioCiio + CoilClll 


= —a 




C100C110 + ClOlClll 


= -/3 




.CiooCioi + CiioCm 





(95) 



with | a |, \/3\, |7| < 2. It is a system of 6 equations in 8 variables. Thus 
the general solutions, for fixed d, live on a 5-dimcnsional manifold. A 
particular 3-dimensional submanifold is obtained by a — (3 = 7 = 0. 
In such a case 

C000C100 = -CooiCioi 
CoioCiio = — C011C111 
C000C010 = —CooiCoii 
C100C110 = —C101C111 
C000C001 = -CoioCoii 
.CiooCioi = — C110C111 

For example, the following MMES is an element of that manifold 
when d = 



1 



= -^=( - |000) + |001) + |010) + |011) 
+ |100) + |101) + |110)-|111)). 



(97) 



As in the case of 2 qubits, this is an example of perfect MMES with 
uniform amplitudes = and real phases Cfe € {—1,1}, with 

kex 3 . 

2. If d = and Cj = c (i e S) with c e [—1/3, 1], one gets 



- (1 + c(cr 1( T2 + 0- 2 0-3 + 0-30-1)) 



Thus, 



|(1 + 3c) for k e {000, 111} 
|(1 — c) otherwise. 

(98) 



+ 3c)(l - c) (C000C100 + C011C111) 

+ (1 - c) (CooiCioi + CoioCiio) = 
v /(l + 3c)(l-c) (C000C010 + C101C111) 

+ (! - c ) (CooiCoii + C100C110) = 
1/(1 + 3c) (1 -c) (C000C001 + C110C111) 

+ (1 - c) (CoioCoii + C100C101) = 



(99) 
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(a) If c = 1 the phases are arbitrary and the perfect MMES is 

1 



V2 



(Cooo|000)+Cin|lll» : 



(100) 



1, one obtains the GHZ 



As a particular case, when Cooo 
state [37]. 

(b) For c < 1, the solutions live on a 5-dimensional submanifold. Note 
that if one tries a solution for which the phases are independent of c 
one gets 

CoooCioo + CoiiCm = 
CoooCoio + CioiCm = 
CoooCooi + CiioCm = 
CooiCioi + CoioCno = 
CooiCon + CiooCno = 
.CoioCon + C100C101 = 



(101) 



that is, 



C100C011 
CooiCno 



; C010C101 = CooiCno 
C010C101 = 
Cooidio + CiooCon = 
k CoioCioi + CiooCon = 

which has no solutions. 



-CoooC 



in 



(102) 



5 Uniform MMES 

According to Theorem [15l a perfect MMES has a population probability 
vector in the computational basis given by (|69p . whose marginals on maxi- 
mal subsystems are all uniform. In particular, a uniform probability vector 
is compatible with a perfect MMES. In this Section we will focus just on 
this class of states, that have uniform amplitudes 

n = W\ = VW) = i/Va, Vk e x n , (103) 

and depend only on N = 2™ phases. 

Definition 14 (Uniform states) A state \tp) € Tis of the form 

l^) = ^E^I fc )' C=&)6T N , N — 2™, (104) 

is said to have uniform amplitudes in the computational basis. A state with 
uniform amplitudes in the computational basis is also called a uniform state. 

First of all, we have a complete characterization of uniform maximizers 
of the potential of multipartite entanglement. 
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Theorem 18 (Fully factorized uniform states) The fully factorized 
states with uniform amplitudes, z — Q/\fN, have = IlieS ' with 

Q i e§\k£X n . 

Proof The result is an immediate consequence of Theorem [5J by observing 
that z k = Uies K ld 4 = IUs K I Il i6 s C4 = Ck/VN. □ 

The various expressions of purity of a bipartition (A, A) considered in 
Section [2] simplify for uniform states. In particular, by plugging (|104() into 
HMD we find 

Theorem 19 (Purity for uniform states) Consider a state with uniform 
amplitudes in the computational basis z — £/\/]V , with £ E T N . Then for 
any bipartition (A, A), 

ta(C) = — — jr^- 1" ^2 ^ ^ Re (Cfe (mi Cfce;©m (k® m ) > 

kexs lex? m£X f 

(105) 

where N A = 2 nA , and N = 2™. 



Proof When z = £/yN, in the first three sums of (f2"6"| all terms are equal 
to 1/N 2 . Their number, according to RemarkCS is Af£t +Mot = N(Na + 
Na — 1), and the result follows. □ 

Remark 25 Note that the first term on the right-hand side corresponds to 
the average entanglement for typical states 22,26,31,32,39], whose phases 
are uniformly distributed on the torus T N . Thus, the combination of phases 
in the second term can increase or reduce the value of the purity with respect 
to the typical one (at a fixed bipartition). 

Finally observe that, by setting ( k = e lVk , with ip k E [0,2%), k E X n , 
one gets 

N A + N A -1 
ta(C) = jj 

+ ftf2 X! X! COS ^ k ~ ^ k ® 1 + Vh®l®m ~ <Pk®m) -(106) 

kex s lex? mex? 

Remark 26 Note that for a uniform fully factorized state, since = Ck A Ck- > 
with Cf 4 - ILe a (l > Vfc € X n , one gets 

(107) 

Vfc e X s , V/ € X A , V/ E X a . Therefore, all terms of the sum in (|TU5J> are 
1, and 

tta(C) = iV - 4+ ^~ 1 + - 1)(N A - 1) = 1, (108) 

as it should. 
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The counterpart of Theorem [19] for the potential of multipartite entangle- 
ment is stated in the following 



Theorem 20 (Potential for uniform states) // the Fourier amplitudes 
in the computational basis are uniform, z — Q/^/N , with Q 6 T N , then the 
potential of multipartite entanglement reads 



tme(C) 



N A + N A - 1 
N 



+ ~fp, E E 5(^ m ;[«/ 2 ]) R- e (CfcCfcei®mCfc®iCfe©m) 



keX n i.tneJf! 



(109) 



where N A = 2^ n / 2 \ N A = 2^ n+1 ^ 2 \ and N = T 



Proof When z = (,/\fN, in the first two sums of |46|) all terms are equal to 
1 /N 2 and one obtains 



^.(O = 4 + 4 ]T 9&\A [n/2] 



N N 



lex: 



72 E E g ( l > m > [ n / 2 D Re Cifceiem Cmi Cfee™) 



N 2 



We get 



E 

l&X" 



g(\l\An A )=J2 E ^(|i|)ff(*.0;n A )= E ' UM;tu) 



lEX? 0<s<n 



and, from lH3l 



;(s,0;n A ) = - 



1 /n 

2 is 



s 



Thus, 



2 ]T ff(|i|,0;n A ) 



E 

Ks<r 



( n A 

V S 



= 2 nA + 2 nA - 2, 



and, by setting n A = [n/2], the result follows. □ 



We will now use Theorem [20] and look for the uniform minimizcrs of the 
potential of multipartite entanglement. 
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5.1 Two qubits 

For two qubits n = 2, we have N = 4, Na = Na — 2 and (|109| becomes 
g ^ 

k£X 2 i.meXJ 

= 4 + g X! ■ 9 ( 1 ' 1; ^ Re ^ fc © 01 C fc ® 10 ) • ( U0 ) 

kex 2 

From (23]) we get §(1, 1; 1), hence 

^me(C) - I + \ Re (CoodiCoiCio) ■ (HI) 
Uniform perfect MMES are solutions of the equation 

7Tme(C) = ^ (£T 4 , (112) 

that is 

CoodiCoiCio - -1, (H3) 

which yields 

life) = 5 (CoolOO) + Coi|01) + Cio|10) - CooCoiCio|H» • (H4) 

In this degenerate case, multipartite entanglement coincides with bipartite 
entanglement, and this state is obviously equivalent, up to local unitaries, 
to a Bell state. A particular subclass is formed by uniform perfect MMES 
(|114[) with real phases £ 6 { — 1,+1} 4 . Their number is 2 3 and has been 
already found by using a probabilistic approach. See (|88[) . 



5.2 Three qubits 

For n — 3 qubits, N = 8, Na = 2, Na — 4, and one must look for the 
solutions of 

7r ME (C) = ~, CeT 8 , (lis) 

where, from (|109p . 

ttme(C) = g + X! X! m ; 1) R- e (Cfe Cfcei®™ Cfce; Cfce™) (H6) 

Due to the constraint Sq(1 A m) in the coupling function g, see Theorem llOl 
one can easily see that the pairs that yield nonvanishing contributions to 
the sum are 

(J,m) = (p(001),p(010)), (I,m) = (p(001),p(110)), P eC 3 , (117) 
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and the pairs obtained by exchanging I and m, where C3 C V3 is the sub- 
group of the 3 cyclic permutations defined in (|57|) . Therefore, 

ttme(C) = g + g2 X! X! [vt 1, 1; 1 ) Re (^ fc Cfeep(on) Cfcep(ooi) Cfe® P (oio) 

+ s(l,2;l) Re(C fe Cfc ep(in) Cfcep(ooi) Cfcep(iioj 



yr^T [2 Rc (Cp(fc) Cp(feeoii) Cp(fceooi) Cp(fc©oio) 



5 

+ Re (Cp(fe) Cp(fe©m) Cp(fe©ooi) Cp(feffiiio) )(H8) 



since from (H31 we get 1; 1) = 1/3 and g(l,2; 1) = 1/6. By performing 
the sum over we finally obtain 



7TMe(C) = g + 4g X! [ 2Re (Cp(000)Cp(011)Cp(001)Cp(010)) 



5 1 

+ 2Re(C p (iii)Cp(ioo)Cp(iio)Cp(ioi)) 



+ R e (Cp(000) Cp(lll) Cp(OOl) Cp(llO) ) 

+ Re(Cp(oio)Cp(ioi)Cp(oii)Cp(ioo)) ■ (119) 

There are 3 families of solutions, living on the following 5-dimensional 
submanifolds 

M p = {(Cfe) £ T 8 I Cp(ooo)Cp(ni)Cp(ooi)Cp(iio) = +1, 

Cp(010)Cp(101)Cp(100)Cp(011) = +1) 

Cp(ooo)Cp(oii)Cp(ooi)Cp(oio) = P ^ C 3 . (120) 

Indeed, if C € M p it is an easy task to see that 

)Cg(000)Cg(lll)Cg(001)Cg(110) = a p- 1 q(l) 
Cg(010)Cg(101)Cg(100)Cg(011) = &p-ig(l) q £ C 3 (121) 

C?(000)Cj(011)Cg(001)Cg(010) = c p- 1 q(l) 
,Cg(lll)C(}(100)Cg(110)Cg(101) = d p -i q (i) 



where 

a = (+1, —a, —a) 
b = (+1, —a, —a) 
c= (-l,+a,-l) 
[d= (-i,+ a ,-l) 



(122) 
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with aG§' arbitrary. Therefore, the sum in 7Tme(C) reads 

Rc (a + b + 2c + 2d) p _ lq{1) = Y Re (a + b + 2c + 2d) i 

p£C 3 i£S 

= 2 - 8 + 4Rca - 4 Re a = -6, (123) 

yielding 7r ME (C) = 1/2. 

Note that, in agreement with Theorem IT2| 7Tme(C) contains Af^ = 12 
distinct terms that depend on phases, 6 of which are double weighted. The 
above solutions force 2 terms to the value + 1, and 4x2 terms to the value 
= — 1. The remaining ones are symmetric around and cancel. 

The corresponding uniform perfect MMES are 

= -^(Cp(ooo)b(OOO)) +C P (ooi)b(001)) +C P (oio)b(010)) 

-Cp(ooo)C P (ooi)Cp(oio) b(on)) + Cp(ioo) b(ioo)) 

-Cp(000)Cp(001)Cp(100) |JP(101)> + Cp(HO) |p(H0)> 
+Cp(ooo)Cp(ooi)Cp(iio)b(lll))), ?eC 3 . (124) 

At present we do not know whether there exist other classes of uniform 
perfect MMES than (|124|) . Numerical evidence seems to corroborate the 
conjecture that (|124j) describe all uniform perfect MMES, but we could not 
prove it. 

5.2.1 Real uniform MMES Let us now look for uniform perfect MMES 
whose phases are all real, i.e. C, € { + 1, l} 8 - A necessary condition is that a 
is real, a G { — 1, +1}- In particular, it is an easy task to prove that a = — 1 
iff a, b, c, d in (|122p are permutation invariant, iff 

a = & = (1, 1, 1), c = d= (-1,-1,-1). (125) 

Thus a = — 1 characterizes the 4-dimensional intersection 

M* = M s o n M s i = M s i n M S 2 = M S 2 n M s o = P| M p . (126) 

P6C3 

On the other hand, a = +1 determines the following three nonintersecting 
4-dimensional submanifolds 

N p = M p n {(Cfc) G T 8 |a = Cp(ooo)Cp(no)Cp(oio)Cp(ioo) = +1}: P € C 3 . 

(127) 

Therefore, all uniform perfect MMES with real C belongs to one of the above 
nonintersecting manifolds, namely 

{real uniform MMES} C |J N p U M*. (128) 

p£C 3 
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Thus the total number of real uniform perfect MMES is 4 x 2 4 = 2 6 . They 
are given by 

M = ^ Ce{-M} 8 (129) 



with 



where 



C000C001C010C011 = x j 
C000C001C100C101 = Vj 
CoooCoioCiooCno = z j 
. CooiCoioCiooCm = w ji 



1 < j < 4, (130) 



z = (-l,-l,-l,+l), y= (-1,-1, +1,-1), 

z =(-!,+!, -1,-1), «;=(-l,+l,+l,+l). (131) 



5.3 n > 3 qubits 

For n = 4 qubits, TV = 16, Na = = 4, and a brute force enumeration 
shows that the minimum value of the potential of multipartite entanglement 
in the class of real uniform states, is 

min {tt^ (0 | C G {-1, 1} 16 } = ^ > \- (132) 
In fact, there are 1056 minimizers, among which, there is, e.g. 

C = (-i,-i,-i,-i,-i,-i,+i,+i,-i,+i,-i,+i,+i,-i,-i,+i). 

. (133) 

There is numerical evidence that 1/3 is the minimum of the multipartite 
entanglement, and thus it is not an artifact of the restriction to real uniform 
states. In fact, it has been proved that for n = 4 the minimum of 7Tme is 
strictly larger than 1/4 [T31l2"B] . but still its value is unknown [5]. This is 
a first example of frustration among the bipartitions, that prevents the 
existence of a perfect MMES: the requirement that purity be minimal for 
all balanced bipartitions generate conflicts already for n — 4 qubits. 

For n = 5 and 6, the expressions become more complicate. Here, we will 
not discuss this cases. We will only exhibit two real uniform perfect MMES, 
solutions to 

4(0 = j Ce{-1,+1} 32 (134) 

and 

4 6 e(0 = ^ Ce{-1,+1} 64 , (135) 

respectively. Therefore, interestingly, frustration is present for n = 4 qubits, 
while it is absent for n — 5 and 6. 
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For example, a 5-qubits real uniform perfect MMES is defined by Eq. 
(|104[1 with the following set of phases 

C = (+i,+i,+i,+i,+i,-i,-i,+i,+i ) -i,-i,+i,+i,+i,+i,+i, 
+1,+1,-1,-1,+1,-1,+1,-1,-1 ) +1,-1,+1,-1,-1,+1,+1) 

(136) 

and can be shown to live on a 7-dimensional manifold, while a 6-qubits real 
uniform perfect MMES has the following set of phases 

— 1, +1, -1, -1, -1, +1, +1, +1,-1, +1,-1,-1, +1,-1, -1,-1, 
+ 1,-1, -1,-1,-1, +1,-1, -1, +1,-1,-1, -1,+1,-1,+1,+1, 
+ 1,+1,+1,-1,+1,+1,-1,+1,-1,-1,-1,+1,+1,+1,-1,+1). 

(137) 

By using the theory of quantum weight enumerators and quantum codes 
[4211430 44]. it has been proved that [45] 

min{ 7 r^(z)|zeS 2 " +1 ^ 1 } > 2-["/ 2 l, for n > 8, (138) 

and thus there is frustration among the bipartitions that prevents the exis- 
tence of a n-qubit perfect MMES, for n > 8. The case n — 7 is still open. 
There is numerical evidence that it is frustrated too, but no conclusive ar- 
guments. 

Summarizing, perfect MMES exist for n = 2, 3, 5, 6 and, possibly, for 
n = 7. For n = 4 and n > 8 there is frustration and the minimum of 
the potential of multipartite entanglement is strictly larger than 2~[™ //2 l. 
Interestingly enough, in the cases considered (n < 6) we have shown that the 
(conjectured) minimum of the potential is attained by uniform states with 
real phases. In such a case, in order to study the structure of multipartite 
entanglement in a quantum state of n qubits, and in particular the minima 
of its potential, one can instead consider the simpler system of classical 
sequences C £ { — 1, +1} 2 " of 2™ bits, with Hamiltonian 7r^g(£). 

6 Conclusions 

In this paper we have studied the properties of the potential of multipartite 
entanglement and of its minimizers, the MMES, for a system of n qubits. 
In particular our focus has been on perfect MMES, that saturate the lower 
bound of the potential, and by using a probabilistic approach, we have 
proven a theorem on the structure of their population probability vectors. 
This allowed us to consider a particular simple class of solutions, those with 
uniform population. We have shown by explicit construction that (apart for 
the case n = 7 which is still open, but probably is frustrated) there always 
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exist uniform perfect MMES with real phases, a class of states that can 
be mapped to the classical binary sequences of length 2™. In fact, we have 
shown that also for n = 4, the lowest number at which frustration occurs 
and hinders the existence of perfect MMES, the (conjectured) minimum 
of the potential of multipartite entanglement is attained by uniform states 
with real phases. This represents a great advantage, because in this situation 
one can investigate the structure of quantum multipartite entanglement by 
studying the simplest problem a classical Hamiltonian defined on binary 
sequences. 
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